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ABSTRACT 

Modeling artificial scanning electron microscope (SEM) and scanning ion microscope images has recently be- 
come important. This is because of the need to provide repeatable images with a priori determined parameters. 
Modeled artificial images are highly useful in the evaluation of new imaging and metrological techniques, like 
image-sharpness calculation, or drift-corrected image composition (DCIC). Originally, the NIST-developed ar- 
tificial image generator was designed only to produce the SEM images of gold-on-carbon resolution sample for 
image-sharpness evaluation. Since then, the new improved version of the software was written in CH — h program- 
ming language and is in the Public Domain. The current version of the software can generate arbitrary samples, 
any drift function, and many other features. This work describes scanning in charged-particle microscopes, which 
is applied both in the artificial image generator and the DCIC technique. As an example, the performance of 
the DCIC technique is demonstrated. 

1. INTRODUCTION 

Computational scanning electron microscopjffl through rapid artificial image modeling is gaining importance. 
It is a useful tool for evaluation of imaging and metrology methods, since real SEMs or other charged-particle 
microscopes cannot always provide repeatable images. For example, it is virtually impossible to obtain two real 
SEM images that only differ in random noise. This is usually caused by many perturbing factors like drift, 
sample charging, or electro-magnetic fields. The artificial image generator is capable of modeling all important 
effects in a deterministic way. One can a priori choose the drift function, the type^ - and magnitude and type 
of noise, the charged-particle-beam profile, etc. That being the case, computer generated artificial images may 
be input to the imaging and metrological techniques and the results compared to the chosen parameters, hence 
indicating the performance of given techniques. None of these is possible with the real images, where these effects 
are present there, but all are random and often even unknown. 

An advanced version of the artificial SEM image generatoi^'^ has been released as a public-domain software. 
It is implemented as a library written in CH — h. This also allows for linking with programs written in many other 
programming languages. The software works in Linux, Mac OSX. Windows, and very probably in other UN*X 
systems as well, however, the latter has not yet been tested. For faster and easier designing of calculations, LucP 
scripting was implemented. Lua is a scripting language originally designed for data-entry applications. These 
days it is mostly employed in computer games. It is one of the simplest and fastest scripting languages available. 
A simple graphical user interface (GUI) has been written mainly for demonstration. One can very easily generate 
images of two types; gold-on-carbon resolution sample and periodic semiconductor cross structures. The GUI 
depends on wxWidget^ library which is multiplatform as well. 

One of the techniques that have been tested with modeled images is the drift-corrected image composition 
(DCIC) which outputs significantly more accurate images than the traditional imaging techniques. This is 
necessary for sub-nanometer-scale metrology, since the conventional "slow-scan" and "fast-scan" techniques pro- 
vide images that are often distorted or blurry. The DCIC works with frames that are taken as quickly as the 
capabilities of the instrument permit. Physical drift causes displacement between each couple of frames. This 
displacement is searched for with cross-correlation. Since the quickly acquired frames are usually extensively 
noisy, a noise reduction is a part of the DCIC technique. 

* Contribution of the National Institute of Standards and Technology; not subject to copyright. Certain commercial equipment 
is identified in this report to adequately describe the experimental procedure. Such identification does not imply recommendation or 
endorsement by the National Institute of Standards and Technology, nor does it imply that the equipment identified is necessarily 
the best available for the purpose. 



2. DRIFT DISTORTION 

In the scanning microscopes, the image is formed by scanning across the sample in a raster pattern. Intensity 
value is acquired at each location on the sample. In digital scanning microscopes, that corresponds with a pixel 
in the image. The intensity value ^(r) depends on the landing position of the electron beam r. Most SEMs use 
the raster pattern. Let the raster pattern be defined by the time-dependent vector function: 
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where t is time, M is a single-pixel step length, x and y are column and row indexes in the SEM image, and 
are the unit vectors in x- and y-direction, tu is the pixel-dwell time, is the dead time between acquisition 
of two pixels, tj is the time needed to move the beam to the beginning of the new line. \_q\ is a symbol for the 
floor(g) function as used in programming languages. X and Y are the pixel-width and pixel-height of the SEM 
image. 

Let the SEM imaging be defined as a relation between the intensity map of the sample ^(f) and the SEM 
image I{x, y): 

i{xit),y{t))=mm)- (4) 

The relation between / and ^ may in practice be very general. For simplicity, lot K he a constant in this 
manuscript, since this does not affect generality of the DCIC technique. In the ideal case: r(t) = T>(t); however, 
drift and space distortions are always present in scanning microscopes and they can significantly affect the 
position r: 

r{t)=frit)+Dd{t) + D,{fr)- (5) 

The space distortion Dg is constant in time and may be compensated for, when its function is known. This dis- 
tortion may be caused by non-linearities in deflection amplifiers and is significant mostly at low magnifications. 
On the other hand, the drift distortion is changing in time, its function is usually unknown, and it may 
extensively affect the high-magnification images. The drift distortion may arise from several sources; e.g. trans- 
lational motion of the sample, tilt or deformation of the electron-optical column, outer forces and vibrations, 
or temperature expansion. High-magnification images are very sensitive to drift distortion, since microscopic 
displacements, tilts, or temperature changes can easily cause nanometer distortions and displacements, which 
can significantly impair the SEM image and its usability for nanometer-scale measurements. 

The drift-distortion function is generally unknown, however, since it characterizes motion of physical bodies, 
it must be continuous and thus square-integrable. Therefore, drift-distortion function may be Fourier-series 
expanded: 
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where c„ are the (complex) Fourier coefficients, U is the overall energy of the drifting system. Since U is limited, 
for high n the coefficients c„ must be nearing zero. In practice, for frequencies higher than 200 Hz, c„ correspond 



2 



Figure 1. Series of artificial images of a semiconductor structure composed using the traditional "fast-scan" technique. 
Compositions of 2, 4, 8, 16, 32, 64, 128, 256, and 512 frames (from the top left). Images are normalized. 

to noise only and are negligible. Therefore, the Dcd{t) can be written: 

N 

Dcd{t)^ J2 Cne-^"*, (9) 

n=-N 

where A'' represents the highest significant angular frequency. 
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Figure 2. Series of artificial images of a semiconductor structure composed using the DCIC method. Compositions of 2, 
4, 8, 16, 32, 64, 128, 256, and 512 frames (from the top left). 



3. "FAST-SCAN" IMAGING 

The imaged intensity signal in the SEM always contains noise. The intensity function is a superposition of a real 
signal and noise: 

ar,t)=Ur)+Ut), (10) 

where is the position-dependent real signal and ^„ is the time-dependent noise. ^„ is a superposition of all noise 
contributions present in the SEM: Poisson noise originating from the electron source and the secondary emission, 
the noise originating from the amplifier and electronics, quantization-error noise, etc. Due to the central limit 
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theorem, it is legitimate to suppose that the mean value of this noise is zero: 



<Ut)>=0. (11) 

In order to obtain a SEM image with a desired level of noise, the overall pixel dwell-time to must be sufficiently 
high. Unfortunately, the electron yield is usually low and the overall pixel-dwell time must often be set to times 
ranging from tens to several hundreds of /is. 

In the SEM, there are two common methods to achieve this, i.e "slow-scan" and "fast scan", while the latter 
is useful for metrological application. 

"Fast-scan" is one of the common imaging methods in SEMs. The image is composed from multiple (N) 
frames, for which averaging is the mostly applied technique. The frames are acquired with the lowest possible 
pixel-dwell time to- The image pixel value is an average of corresponding frame- pixel values: 

Ik{x{to),y{to)) = KUr{to + ktf)) + 

+ K^nito + ktf), (12) 
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tf = Y{Xtp + t,) + tjj, (14) 

tf is a, time period between beginnings of acquisition of two following frames, tjj is the dead time between the 
end of acquisition of one frame and beginning of the next one. Considering Eq ()11|) . the higher the lower 
noise level is present in the image. The required noise-level thus determines the number of composed images N. 
For high N: 

N-l 

J2^n{to+ktf)^0. (15) 
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Because the scanning raster pattern is constant for all frames, 

fr{to + ktf) ^f{to). (16) 

Eq (fT2|) may be expanded: 

I{x{to),y{to)) = ^E6[n.(io) + 
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+ Ds{rrito)) + Dd{to + ktf)]. (17) 

With current SEMs, the frame-acquisition time tf can be much lower than the period of even the highest drift- 
distortion frequencies. The drift-distortion within the single-frame acquisition time is then minimal. However, it 
becomes significant during acquisition of the whole image, especially, when the dead times tjj are prohibitively 
high, which is the case even with some current instruments. 

4. DRIFT-CORRECTED IMAGE COMPOSITION (DCIC) 

The "fast-scan" method may be significantly improved using drift-distortion correction, when the images are 
acquired quickly enough. Since the space-distortion Ds is much less pronounced and much smaller that the 
drift-distortion Dd at very high magnifications, it will be neglected from now on. The Eq ([T7|) then becomes: 

H^,y) - J7 ^^^^r{r) + Dd,], (18) 
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Figure 3. Error distribution of the displacement detection. Artificial SEM image of a periodic semiconductor sample was 
used. 

The image is in this case the mean value of N displaced images. 

Under certain conditions, it is possible to find the displacement vectors of the images, which are equal to 
the drift-distortion values Ddk- The drift-distortion then may be compensated for, which allows for acquisition 
of a corrected, more accurate image. One possible approach is a cross-correlation-based displacement detection, 
which is used in the DCIC technique. The maximum of the cross-correlation function is searched for. Its position 
is equal to the searched displacement vector D^k- 

In the DCIC technique, the cross-correlation with noise reduction is applied. This is necessary, because 
the quickest-acquired images are usually very noisy and the peak in the cross-correlation function becomes 
overridden by numerous other peaks, corresponding to random correlation of noise. This often makes finding 
the displacement vector impossible. This issue can be tackled by low-pass frequency filtering performed in the 
frequency domain. The cut-off frequency is determined by the filter-radius R. 

Plain maximum search in a discrete function limits the accuracy to a minimum of one pixel. However, in 
the DCIC, the detection of the displacement vector Ddk is performed with sub-pixel resolution. The peak in 
the two-dimensional cross-correlation function is interpolated with a polynomial third-order two-dimensional 
polynomial function and the algorithm then searches for its maximum. 

The technique is very powerful, since it can correct for the drift-related distortions and blur in extremely 
noisy images. (See Figs [T] and [5]) 

5. ACCURACY OF THE DCIC TECHNIQUE 

The accuracy of the detected displacement vector D^k characterizes the accuracy of the DCIC imaging technique. 
Errors in the displacement vector can cause blur. Such blur can under certain circumstances be larger than with 
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application of the original "fast-scan" technique {D^k = 0). In metrological applications, where dimensions are 
measured from the images, the drift-related displacement is the main source of errors. 

The artificial SEM images have been successfuly used to evaluate accuracy of the DCIC technique. The 
artificial-image generator is, unlike any other source of SEM images, capable of modeling all necessary charac- 
teristics for this application, e.g. arbitrary drift functions, dead times, arbitrary types of samples, etc. 

A performance characteristics must be chosen to investigate the limits of an imaging technique. Application 
of the standard deviation of the displacement vector would be a good candidate, if the distribution of errors was 
Gaussian. In order to find this out, a large set of artificial images (500 000) randomly differing in displacement 
and noise has been applied to find the error distribution of the displacement detection. The DCIC technique 
has processed all generated frames and has output corresponding displacement values. The two-dimensional 
histogram of these values forms the resulting distribution, which is shown in the Fig [3] These data have clearly 
indicated that the error distribution is not (always) Gaussian. Using standard (Gaussian) error processing has 
therefore been unsuitable and thus we have chosen the mean error 6o as the performance characteristics. 

N-l 

fc=i 

where Sok is the error of the displacement vector Ddk and N is the number of frames. Since the correct 
displacement vector Dd is known (it is determined by the artificial-image generator), 

Sok = \Ddk - D,k\. (21) 

The performance of the DCIC technique is obviously limited, because noise, blur, contrast, and other param- 
eters affect it significantly. For instance, if the frames were extremely blurred and the cross-correlation maximum 
would be overly wide and the mean error of the displacement vector would be excessively high. It is therefore 
useful to find the dependences of 6d on noise and blur and provide a set of limiting parameters. 

The dependence of the mean error of the detected displacement 6d on noise and blur have been both inves- 
tigated with application of artificial images. Gaussian noise and Gaussian blur have been chosen for simplicity, 
although the type of noise and the blur profile may be arbitrary. For every step in noise and blur, 5000 artificial 
images of the gold-on-carbon resolution sample have been generated and processed by the DCIC algorithm. The 
results of these tests are shown in Figs 0] and [HI For reference images showing different magnitudes of Gaussian 
noise see Fig [51 These tests demonstrate the capability of the DCIC technique to find the displacements with 
sub-pixel accuracy. In the noise test, this is maintained up to the ag ~ 8, which roughly corresponds to signal 
to noise ratio (SNR) around 0.1 and the dependence is almost linear. The dependence on blur indicates that the 
sub-pixel accuracy is sustained up to cTf, = 14. 

6. CONCLUSION 

Modeled artificial SEM images were first employed in assessment of the image-sharpness calculation technique^ 
and have been adopted as a part of the developed international standard for image sharpness. Since then, a 
new highly improved version of the software was written. This version supports arbitrary non-overlapping two- 
dimensional samples, rigorous generation of Poisson and Gaussian noise, arbitrary drift functions, dead times and 
other features. Scripting in Lua scripting language was implemented to make the calculations easier to design. 
This new tool was then used in evaluation of the new imaging technique of DCIC. By finding dependence of the 
error in detection of the displacement on noise and blur, the sub-pixel accuracy was demonstrated even for high 
magnitudes of noise or blur. This makes the DCIC and modeling of microscope images useful and important 
tools for nanoscale metrology and nanotechnology. 
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Figure 4. Evaluation of the DCIC technique. Dependence of the mean displacement detection error on the magnitude (ug) 
of Gaussian noise. Each represents 5000 artificial images of the gold-on-carbon resolution sample sized 512x512 pixels. 
The error-bars denote the standard deviation of the displacement vector detection error. 
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Figure 5. Gaussian-noise scale. Artificial images of the gold-on-carbon resolution sample with Gaussian noise of different 
magnitudes. 
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Figure 6. Evaluation of the DCIC technique. Dependence of the mean displacement detection error on Gaussian blur 
(db). This blur simulates the effect of the charged-particle-beam profile. Each represents 5000 artificial images of the 
gold-on-carbon resolution sample sized 5f2x512 pixels. The error-bars denote the standard deviation of the displacement 
vector detection error. 
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